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Abstract 



Background: Analysis of data involving nuclei far from stability often requires optical potential 
(OP) for neutron scattering. Since neutron data is seldom available, while proton scattering data 
is more abundant, it is useful to have estimates of the difference of the neutron and proton optical 
potentials. This information is contained in the isospin dependence of the nucleon OP. Here we 
attempt to provide it for the nucleon-^^'^Pb system. 

Purpose: The goal of this paper is to obtain accurate n+^'^^Pb scattering data, and use it, together 
with existing p+^^^Ph and ^^^Ph{p, n)^'^^Bij^g data, to obtain an accurate estimate of the isospin 
dependence of the nucleon OP at energies in the 30-60 MeV range. 

Method: Cross sections for n+^'^^Pb scattering were measured at 30.4 and 40.0 MeV, with a 
typical relative (normalization) accuracy of 2-4% (3%). An angular range of 15 to 130 degrees 
was covered using the beam-swinger time of flight system at Michigan State University. These 
data were analyzed by a consistent optical model study of the neutron data and of elastic p+^^^Pb 
scattering at 45 MeV and 54 MeV. These results were combined with a coupled-channel analysis 
of the ^'^®Pb(p, n) reaction at 45 MeV, exciting the 0'^ isobaric analog state in ^"^^Bi. 

Results: The new data and analysis give an accurate estimate the isospin impurity of the nucleon- 
208p]-, Qp 3Q_^ MeV, caused by the Coulomb correction to the proton OP. The corrections to 
the real proton OP given by the CH89 global systematics was found to be only few percent, while 
for the imaginary potential it was over 20% at the nuclear surface. Based on the analysis of the 
measured elastic n+^^^Pb data at 40 MeV, a Coulomb correction of similar strength and shape 
was also predicted for the p-|-^*^'^Pb OP at energy around 54 MeV. 

Conclusions: Accurate neutron scattering data can be used in combination with proton scattering 
data and (p, n) charge exchange data leading to the IAS to obtain reliable estimates of the isospin 
impurity of the nucleon OP. 

PACS numbers: 25.40.Dn; 24.10.Ht; 24.10.Eq 
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I. INTRODUCTION 



At energies below 100 MeV, the attraction between a neutron and a proton is stronger 
than that between two protons or two neutrons. Consequently, the average interaction of 
a proton with an > Z nucleus is stronger than that of a neutron. In other words, the 
nuclear interaction between an incident nucleon and a target with non-zero isospin has an 



isospin dependent part. For the nuclear part of t 
the isospin dependent term is, in the Lane form 



le nucleon- nucleus optical potential (OP), 



t.T 



(1) 



where t is the isospin of the incident nucleon and T is that of the target A. The second 
term of Eq. ([T]), known as the Lane potential, contributes to both the elastic {p,p) and 
{n,n) scattering as well as to the charge exchange (p, n) reaction {2!, 3|. A knowledge of Ui 
is of fundamental interest for studies of nuclear phenomena in which neutrons and protons 
participate differently (isovector modes). Many previous estimates of f/i (see, for example, 
Refs. |#^6|) involved a comparison of the nucleon OP from a range of nuclei with different 
values of the asymmetry parameter e = {N — Z)/A ot the Distorted Wave Born Approxima- 
tion (DWBA) analysis of {p,n) reactions exciting the isobaric analog state (IAS). However, 
these approaches are subject to serious uncertainties. For example, in the comparison of 
elastic nucleon scattering from different nuclei one must make assumptions j4| about the 
variation of nuclear geometry with A and e. Moreover, the contribution of the Lane poten- 
tial f/i to the elastic nucleon scattering cross section is relatively small [3, 8|. In contrast, 
Ui entirely determines the (Fermi-type) AJ^ = O"*" transition strength of the {p, n) reaction 
exciting the IAS; this reaction is, therefore, a sensitive probe of the isospin dependence of 
the proton-nucleus OP jsl. However, for {p,n) reactions, a change in Re Ui can be ap- 



10|, and the determination of f/i remains 



proximately compensated by a change in Im Ui 
ambiguous without additional information or constraints. 

It is in principle possible to avoid these uncertainties by extracting Ui from a consistent 
study of the elastic proton and neutron scattering and the charge exchange {p, n) reaction 
on the same target nucleus, at the same energy. We recall here briefly the consistent isospin 
coupling scheme jsj] for the elastic nucleon- nucleus scattering and charge exchange (p, n) 
reaction exciting IAS. For the isospin projections = (A^ — Z)/2 of the target nucleus A 
and Tj = T2 — 1 of the isobaric analog nucleus A, and denoting the neutron- and proton 



scattering states by \nA) and \pA), respectively, the neutron and proton optical potentials 
are given by the diagonal matrix elements of potential ([1]) 

Up = {pA\Un\pA) = Uo- jT^Ui = Uo- eUi, (2) 

Un = {nAlU^lriA) = Uo + jT^i = Uo + eU,, (3) 

Similarly, the transition matrix element or {p, n) form factor (FF) for the charge exchange 
n)y4iAs reaction exciting IAS is 

Fp„ = (ni|4f/i^|M) = |v^f^i = t^i- (4) 

If the neutron and proton optical potentials at a given energy are well determined from the 
optical model (OM) analysis of the corresponding elastic data, then the isovector term of 
the nucleon OP can be obtained directly from Eqs. (EJ-Q as 

U. - (5) 

2e 

Unfortunately, isospin is not a good quantum number in the Coulomb field of the nucleus 
since this field slows the incident proton and affects the strength and shape of the proton- 
nucleus OP. It is necessary, therefore, to add Coulomb corrections AEc to the incident 
proton energy and AUc to Up to separate the main effects of the Coulomb field so that 
isospin is a good quantum number for the remainder of the nucleon OP, namely, 

Up = Uo- eUi + At/c, (6) 
Un = Uo + eU,. (7) 

Then the Coulomb correction term must be determined from 

AUc = Up-Ur, + 2eU, ^ U, - ^ = ^^V^- (») 

In general, Eq. ([8]) has 4 unknowns (the real and imaginary parts of AUc and Ui), that, as 
we shall see, cannot be determined unambiguously even if the complex neutron and proton 
optical potentials are well determined. 

There are "global" systematics of the OP parameters deduced from the extensive OM 



analyses of nucleon elastic scattering, for example, those by Becchetti and Greenlees by 



Koning and Delaroche 



12|, and by Varner et al. 13|. In the work described here, we rely 



on the CH89 global model by Varner et al. The CH89 optical potentials cover a wide range 



of energies and target masses, and are parametrized using Woods-Saxon (WS) forms. The 
resulting systematics are often used to predict the nucleon OP when elastic scattering data 
are not available or cannot be measured, as is the case for many unstable nuclei near the 
dripline. Given the large isospins of dripline nuclei, it is important to estimate accurately 
the isospin dependence of the nucleon OP (or equivalently, the Coulomb corrections to that 
OP) before applying it in studies of nuclear reactions or of astrophysical phenomena. So 
far, the empirical isospin dependence of the nucleon OP has been deduced |lll-ll3l] based 
mainly on the OM analyses of the proton and neutron elastic scattering, adopting some 
simple treatment of the Coulomb correction terms AEq and AUc- 

Before the present measurement of neutron scattering on ^^'^Pb, detailed elastic n+'^^^Ph 
scattering data were available only at energies up to 26 MeV. This energy range, however, 
does not overlap that of much of the precise proton data and furthermore is too small to 
establish clearly the energy dependence of the neutron OP. The measurements at 30.4 and 
40.0 MeV described here have greatly expanded the energy range for the neutron data and 
provided the most accurate and detailed data for neutron scattering above 26 MeV from 
any N ^ Z nucleus. After our measurement, several experiments on the elastic n+^'^^Pb 
scattering have been carried out and the neutron scattering data were measured at 65 MeV 
jl^ . in the forward angular region but over a wide range of the neutron incident energy 
from 65 to 225 MeV [ISj, and at the neutron energy of 96 MeV 16|. Together with the 
present data, one has now a good database of the elastic n+'^^^Ph scattering data to study 
the energy dependence of the neutron OP. In the present paper we show that the accurate 
neutron scattering data can also be used in a consistent analysis of the elastic neutron and 
proton scattering from ^''^Pb and the charge exchange ^''^Pb(p, n)^'^^Bij^g reaction to study 
the isovector part of the nucleon-nucleus OP, and to estimate the Coulomb correction term 
dH]) to the proton OP. 



II. EXPERIMENTAL METHOD 



he measurements were performed using the MSU beam swinger time-of-fiight system 



17H19| as modified for neutron scattering. Neutrons produced by the ^Li(p, n)^Be(g.s. + 
0.429 MeV) reaction scatter from a cylindrical 200.64 gram target (2.40 cm diameter x 
3.90 cm long) of isotopically enriched ^osp]^ (98.69%) |20| and are detected in a liquid 
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scintillation counter with an overall time resolution of about 1.0 nsec. This corresponds 
to an energy resolution for the elastic peak of better than 1.1 MeV FWHM, sufficient to 
resolve the first excited state of ^'^^Pb at 2.61 MeV. Pulse shape discrimination is utilized 
to eliminate the 7-ray background. The neutron detectors are situated in a room separated 
from the swinger vault by a 1.8 m thick concrete wall, except for a hole to transmit the 
target scattered neutrons. Additional shielding against neutrons coming directly from the 
production reaction is provided by a movable 1.1 m long iron shadow bar. A monitor time- 
of-fiight detector is mounted rigidly to the beam swinger so as to measure neutron fiux from 
the production reaction at a fixed angle near 22°. This monitor is used to normalize the 
fiux from run to run. Air scattering background is accounted for by measuring target-in and 
target-out spectra at each angle; a small correction is made to account for the fact that some 
of those air-scattered neutrons originating behind the sample are absorbed by the sample 
on their way to the detector 19|. Relative uncertainties are typically 2-4% but reach 8% 
at a few angles. Observation of the '^Li(p, n) fiux at 0° measures the product of incident 
neutron fiux and detector efficiency, and yields the absolute normalization to within 3%. 
Corrections are made for the dead time, source anisotropy and background attenuation due 



to the sample. Further details of the experimental procedures can be found in Refs. 



18 



19|. 



III. DATA ANALYSIS AND DISCUSSION 



A relatively large target was used to obtain adequate statistical accuracy, making it nec- 
essary to correct the experimental data for the effects of multiple scattering, angle averaging 
and attenuation. Because the cross sections varied rapidly with angle we were concerned 
that the deconvolution procedures generally employed would lead to ambiguities and un- 
acceptably large uncertainties. We elected to avoid these uncertainties by convoluting the 
results of optical model predictions, a straightforward procedure, before comparing them 
with the data in a search routine. For this purpose, Kinney's finite geometry code 2l| was 
incorporated as a subroutine of the optical model search code GIBELUMP j22|. The spin- 



orbit part of the neutron OP was fixed at that used in Ref. [23| for lower energy neutron 
data, and the "average" geometry of the WS potential was taken from Ref. 241 ]. 

During the search, the smeared optical model cross sections were compared with the 



experimental data and the optical model parameters were adjusted until a good fit to the 



experimental data was obtained, thereby fixing the optical model parameters. The cross 
sections and their uncertainties resulting from this procedure were then deduced. These 
cross sections are corrected for multiple scattering, angle averaging, and attenuation, and 
are shown in Fig. [H Fig. H] and Table HI The uncertainties resulting from the finite geometry 
corrections varied with angle, and have a maximum of about 8% at the first diffraction 
minimum. These uncertainties are included in the tabulated cross sections. Because this 
is an unconventional procedure its results were compared with the normal deconvolution 
procedure used at lower energies for the case of Fe at 26 MeV; the two procedures yielded 
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19|. 



close agreement. The analysis process is described in much more detail in Refs. 

of the elastic n+^'^^Pb scattering at 30.4 and 40 



The deduced cm. cross sections 



MeV have been studied in several OM analyses, including the extensive searches for global 
parameters of the nucleon OP 12|. In the present work, a detailed OM analysis of the 
measured elastic n+^^^Pb scattering data was made using the CH89 geometric parameters 
of the nucleon OP 13|. The OM analysis and coupled- channel calculation of (p, n) reaction 



were made using the code ECIS97 written by Raynal [26 1 



A. Coulomb correction to the proton incident energy 

The Coulomb correction to the proton incident energy arises because the proton slows 
down in the repulsive Coulomb field of the nucleus and, hence (because the real part of the 
OP decreases with increasing energy), the real part of Up is more attractive compared to 
that of the neutron OP at the same bombarding energy, even when Ui = 0. Estimation of 
the Coulomb correction of the nucleon OP requires that one first determine the difference in 
the effective proton and neutron incident energies, so that the same isoscalar and isovector 
potentials ?7o(i) can be used to generate the proton and neutron OPs using Eq. ([2]). 

The difference in the proton and neutron energies has usually been assumed either to be 
constant at the average Coulomb energy of the incident proton, AEc = GZe^ /{5Rc) ~ 19 
MeV for a ^o^Pb target Q, or to be energy dependent [12]. In the present work we have 
chosen the CH89 global parametrization for the nucleon OP, which has a simple functional 
form and is quite reliable for the nucleon elastic scattering from medium-mass nuclei at 
energies of 10 to 65 MeV, as the starting point of our OM analysis. The CH89 complex 
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FIG. 1. (Color online) Elastic n+^^^Pb scattering data at 30.4 and 40 MeV, corrected for finite 
geometry effects as described in the text, are compared with the OM fit given by the modified 
CH89 optical potential (see OP parameters in Table ln|) . and represent our best fit to the neutron 
data. 

nucleoli OP is determined as the following energy dependent functional 

U{E) = Uo{E - AEc) T eUi{E - AEc), (9) 

with - and + sign pertaining to the incident proton and neutron, respectively. Thus, the 
CH89 proton and neutron optical potentials determined at energies Ep and En = Ep — AEc, 
respectively, should be fully consistent with the Lane formalism (l2])-([3]). We discuss below 
the extent to which this assumption accurately describes the experimental data. 



8 



TABLE I. Measured cross sections of 2°^Pb(n, n)^°*Pb at 30.4 and 40 MeV. The cross sections 
have been corrected for finite geometry effects as described in the text. 
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To obtain an estimate of AEc based on experiment, we have determined the proton 
bombarding energy at which the slowed proton and the neutron have the same average 
momentum, and thereby the diffraction maxima and minima of the proton and neutron 
angular distributions fall at about the same (average) angles in the forward region {2^ . For 
this purpose, the diffractive oscillation of the elastic n+'^^^Ph data at 30.4 MeV has been 
compared with those of the elastic p+'^^^Ph scattering at 45, 47.3 and 49.4 MeV (see Fig. |2]). 
It is clear that the oscillation pattern of the 30.4 MeV neutron data does not agree with 
that of the 49 MeV proton data, as would be expected from the average Coulomb energy 
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FIG. 2. (Color online) Diffractive oscillation of the elastic n+ Pb scattering data at 30.4 MeV 

25|. 



in comparison with those of elastic p+ Fh scattering data at 45, 47.3 



2|] and 49.4 MeV 



prescription of the CH89 model. The elastic p+^°^Pb scattering data at 45 MeV 2J] do 
have an oscillation pattern at forward angles very similar to that of the 30.4 MeV neutron 
data. This implies that AEc ~ 14.6 MeV, significantly smaller than the value of 19 MeV 
used in the CH89 parametrization. 

Our procedure presumably equalizes the average momenta or wavelengths of the neutrons 
and protons in the regions dominating the scattering. For illustration, we have plotted in 
Fig. |3] the local nucleon (relative motion) momentum K{R) determined from the real folded 
nucleon OP [3] as 



K\E,R) = ^[^c.m. - V{R) - VciR)], 



(10) 



where n is the nucleon reduced mass, V{R) and Vc(-R) are, respectively, the real central 
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FIG. 3. (Color online) Local nuc 
given by the folding calculation 



eon momentum K{R) determined from Eq. (jlOp with the real OP 



13, 



using the density dependent CDM3Y6 interaction 



nuclear and Coulomb parts of the OP (Vc = for the neutron-nucleus system). It can be 
seen from Fig. [3] that the local momentum of 30 MeV neutrons is equal that of 45 MeV 
protons at i? ~ 6 fm, in the surface region of the ^osp]-, target. Because the diffractive 
scattering is dominantly determined by the strength and shape of the OP at the surface, 
the results plotted in Fig. [3] support a Coulomb correction of ~ 14.6 MeV at a proton 
incident energy of 45 MeV. 

Finally, in Fig. H] we show OM calculations with several related optical potentials to 
assess differences in the cross sections they predict. First, the two data sets were compared 



with calculations using the original CH89 parameters [13|, but with Ai^c* set to 14.6 MeV 
at 45 MeV proton energy. Then the potential depths were adjusted, yielding mCH89 OP, 
to give the best fit to the data. For a comparison, OM results given by the complex 



folded potential |9[ calculated with the CDM3Y6 interaction 
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FIG. 4. (Color online) OM description of the elastic n+^^^Pb data at 30.4 MeV and p+^^^Pb data 



at 45 MeV 



2|] given by the CH89 (dotted hne), mCH89 (sohd line), and the folded (dash-dotted 



line) optical potentials. 



imaginary folded potentials were adjusted to the best fit to the data) are shown. 



B. ^"^Ph{p, n)^"^BiJ^g data and the isovector term of the OP 

In the two-channel approximation for the charge exchange (p, n) reaction exciting IAS, 
the total wave function can be written as 



^ = \pA)xpAiR) + \nA)XnA{R), 



(11) 



where the x{R) describe the relative nucleon-nucleus motion. The elastic {p,p) scattering 
and charge exchange A{p,n)AiAs cross sections are then obtained from the solutions of the 
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FIG. 5. (Color online) CC description of the charge exchange Pb(p, n) Bij^g data measured 



at 45 MeV 



29( 1 given by Fpn built upon the isovector part Ui of the CH89 nucleon OP The 



two different curves were obtained with two choices of the 45 MeV proton OP (see Table |II|) 
following coupled-channel (CC) equations jsj 



(12) 
(13) 



Here -ft'p(n) and -E'p(n) are the kinetic-energy operators and cm. energies of the p+A and n+A 
partitions. The proton OP in the entrance {p + A ) channel was determined from the best 
OM fit to the elastic p+^^^Ph scattering data at 45 MeV (see Fig. [Hand potential parameters 
in Table II), while the neutron OP in the outgoing {n + A ) channel was constructed from 
the isoscalar and isovector parts of the nucleon OP using the standard isospin coupling 
scheme js, 9|. Since the energies of isobar analog states are separated approximately by the 
Coulomb displacement energy, the (p, n) transition between them has a nonzero Q value. To 
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account for this effect, the isoscalar Uo and isovector Ui potentials used to construct Fpn{R) 
and Un{R) were evaluated from the CH89 systematics at an effective incident energy of 
E = i^iab — <5/2, midway between the energies of the incident proton and emergent neutron 
[2I. Given the elastic p+^°^Pb scattering data {2^ and charge exchange ^'^^Pb(j), ?7,)^°^Bij^g 
data {29] (both measured at 45 MeV), we were able to test the isovector term f/i of the proton 
OP by comparing the results of the CC calculation for the (p, n) cross section with the data. 
One can see in Fig. |5]that the isovector part of the CH89 nucleon OP 13| accounts very well 
for the {p, n) data, especially when the WS strengths of the proton OP are optimized by 
the best fit to the elastic {p,p) data at 45 MeV (the mCH89 potential in Table HB . This 
result shows that the complex isovector potential Ui given by the CH89 parametrization for 
the proton OP at 45 MeV can be used to estimate the Coulomb correction using Eq. 
and the corresponding neutron optical potential [/„. 



C. Coulomb correction to the proton OP 



After AEc has been fixed to give the same diffraction patterns at the forward angles 
for the proton and neutron elastic cross sections at Ep and En = Ep — AEc, respectively, 
one might naively expect from Eq. (Q that the corresponding proton and neutron optical 
potentials are fully Lane consistent. But this is only true if the only physics involved is that 
due to the energy shift. Yet we know there are other effects, mostly affecting the imaginary 
potential, for example. Coulomb excitation, different Q values for (p, n) and {n,p) reactions, 
and different level structures. In the following we obtain an estimate of the importance of 
such phenomena. 

Given that the isovector potential Ui of the CH89 potential was shown above to give 
a realistic description of the {p,n) data, it is reasonable to use Ui given by the CH89 
systematics to estimate AUc, based on the OM analysis of the elastic p+^^^Pb and ra+^^^Pb 
scattering data measured at 45 and 30.4 MeV. In our OM analysis we used the same WS 



functional form for the nucleon OP as that used by the CH89 systematics 
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13|, so the real 



and imaginary parts of the nucleon OP are determined as 

V{R) = -VJUR,Rv,a,) + ^^fUR,Rso,aso){l.(T), (14) 

W{R) = -W^Us{R, Rw, a^) + 4a^W,^/ws(i?, Rw, (15) 
where f^s{R,Rx,a^) = 1 / {I + exp[{R - R^) / a^]} ■ 

The real part of the proton OP includes the Coulomb potential Vc{R) taken from the CH89 
systematics [isl. The CH89 parametrization gives the isovector part Ui of the nucleon OP 

as 

Ui{R) = Vif„s{R,Rv,av) - i^ayjWi-^f^^{R,Ryj,aw). (16) 

For the proton OP at 45 MeV, the complex strength of the isovector potential is readily 
obtained from the CH89 parametrization as eVi = 2.77 ± 0.17 MeV and eWi = 2.05 ± 0.13 
MeV. The quoted uncertainties were determined from the systematic uncertainties of the 
CH89 global parameters 131 ]. 

In the present OM analysis, we first generated the complex p+^°^Pb and n+^'^^Pb optical 
potentials using the CH89 parametrization but with a corrected AEc = 14.6 MeV for the 
proton OP (see CH89 parameters in Table [TTl)- As a result, the CH89 optical potentials 
for 45 MeV proton and 30.4 MeV neutron are fully Lane consistent, with the Coulomb 
correction taken into account only by using the energy shift AEc in Eq- (E]), and its remnant, 
determined by Eqs. f|T7j) - f|T8l) . is exactly zero. The OM description of the 30 MeV neutron- 
and 45 MeV proton elastic scattering by the Lane consistent CH89 optical potentials is 
shown in Fig. H] as dotted line. One can see that the Lane consistent CH89 OP describes 
the data fairly well, excepting at the backward angles where it fails to follow the oscillation 
pattern of the measured {p, p) data. A much improved OM description of the data has been 
achieved by adjusting the depths of the OP while keeping the same WS geometry as that of 
the CH89 potential (see mCH89 parameters in Table [TTl and solid lines in Fig. |4]). Such an 
adjustment procedure naturally gave rise to a non-zero remnant of the Coulomb correction 
AUc = AVc + iAWc- Assuming the same Ui for the mCH89 OP as that of the original 
CH89 OP, the complex remnant of the Coulomb correction can be explicitly determined 
from Eqs. p.CT. and (fT5l) as 

AVciR) = -AVJM Rv. a.), (17) 
AWc{R) = -AWJUR. Rv^. + ^a^AWs^fUR. R^, (18) 
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TABLE II. OP parameters ([Il])-([l5]) used in the OM analysis of the elastic p+^^^Pb data and 
^_l_208pj-j (jg^^g^^ The radii and difFuseness parameters were given by the CII89 parametrization (in 
fm): Ri, = 7.18, Rw = 7.46, = = 0.69, i?so = 6.73, Oso = 0.63. The incident energy 
E, Coulomb corrections (jl7p -(|18p and potential depths are given in MeV. The real strength of 
the isovector potential was taken as eVi = 2.77 it 0.17 MeV, its imaginary strength eWi and the 
Coulomb correction strengths AV^, AW^, and AWg are given with uncertainties (in brackets) 
estimated from those of our OM fit and the CH89 potential parameters [l^ . 



System 


E 


Potential 


K 


AK 




AW, 




AWs 


eWi 




p+208pb 


54.2 


CH89 


43.71 


0.00 


4.50 


0.00 


6.53 


0.00 


-1.80 (0.11) 


5.90 






mCH89 


40.61 


-1.90 (0.56) 


4.67 


-0.20 (0.30) 


5.03 


-1.88 (0.51) 


-1.80 (0.11) 


5.90 


n+208pb 


40.0 


CH89 


38.17 




4.50 




2.93 




1.80 (0.11) 


5.90 






mCH89 


36.97 




4.87 




3.31 




1.80 (0.11) 


6.11 


p+208pb 


45.0 


CH89 


46.58 


0.00 


3.34 


0.00 


7.43 


0.00 


-2.05 (0.13) 


5.90 






mCH89 


43.27 


-1.94 (0.54) 


3.69 


-0.75 (0.26) 


5.72 


-1.18 (0.52) 


-2.05 (0.13) 


4.22 


n+208pb 


30.4 


CH89 


41.04 




3.34 




3.33 




2.05 (0.13) 


5.90 






mCH89 


39.67 




4.44 




2.80 




2.05 (0.13) 


6.47 



The deduced WS strengths f|T8|) of AUc, further referred to as the Coulomb correction to the 
mCH89 proton OP, are given in Table [ITl with the uncertainties estimated consistently from 
the standard errors of the OM fits and those of the CH89 potential parameters. Although 
the strengths of the spin-orbit potential Vso of the best fit mCH89 optical potentials also 
differ from the original CH89 values, we did not assign this difference to the isospin impurity 
of the OP caused by the Coulomb correction. 

From the strengths of the Coulomb correction AUc to the best-fit mCH89 optical poten- 
tials for the elastic neutron and proton scattering at 30.4 MeV and 45 MeV, respectively, 
we find following ratios of the Coulomb correction to the strength of the OP: AV^/V, ~ 
4.5%, AWJW, ^ 20.3%, and AWs/Ws ^ 20.6%. It is obvious that, within the Lane for- 
malism ([T]), these ratios give us a realistic estimate of the isospin impurity of the mCH89 
nucleon OP. The radial dependence of the Coulomb correction to the real and imaginary 
parts of the central OP is plotted in Figs. [6] and [TJ respectively, where one can see clearly a 
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FIG. 6. (Color online) Real central OP of 45 MeV protons and the corresponding Coulomb cor- 
rection AVc determined from Eq. ()17|) using the mCH89 optical potentials for 30 MeV neutrons 
and 45 MeV protons. The uncertainties were estimated from that give n by the CH89 systematics 



for eAVi and standard errors of the OM fit with the code ECIS97 



261]. 



repulsive character of the Coulomb correction to both the real and imaginary OP. In terms 
of the isospin impurity, we conclude that the Lane formulation ([T]) is accurate to within 
about 4-5% for the real central OP. However, the isospin impurity becomes much larger 
(above 20%) for the imaginary part of the CH89 OP and reaches its peak at the surface. 

In a study of the elastic n+^'^^Pb scattering at 40 MeV similar to that described above 
for 30.4 MeV neutrons, and using the complex CH89 OP Q to predict the elastic p+^°*Pb 
scattering at the higher energies, we find that the elastic proton scattering at 54.2 MeV has 
about the same forward angle diffraction pattern as the measured elastic ra+^^^Pb scattering 
at 40 MeV, yielding AEc ~ 14.2 MeV for the 54 MeV proton potential. There are no 
elastic proton data available near 54 MeV, so we compare in Fig. [8] only the OM prediction 
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for the elastic p+^^^Pb scattering at 54.2 MeV with the elastic n+^°^Pb data at 40 MeV. 
In order to roughly estimate the Coulomb correction for 54.2 MeV protons we have scaled 
the real and imaginary WS strengths of the CH89 OP by the same factors deduced from 
the corresponding strengths of the mCH89 OP compared to those of the CH89 OP for 45 
MeV protons. The parameters obtained are given in Table [TTl as mCH89. We note that the 
imaginary parts of the CH89 and mCH89 proton OPs at 54.2 MeV give total reaction cross 
sections (Jr ~ 2121 and 2050 mb, respectively, and that the latter value agrees quite well 
with the experimental trend: measured values of ctr lie around 2000 mb at proton energies of 



40 to 60 MeV |30[. We have further assumed the same Ui for the mCH89 proton OP at 54.2 
MeV as that given by the original CH89 parameters for 54.2 MeV protons, and deduced the 
Coulomb correction Af/c" to the mCH89 proton OP in the same way as done above for 45 
MeV protons. We find a behaviour of AUc very similar to that found for 45 MeV protons, 
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FIG. 8. (Color online) OM description of the elastic n+^o^Pb data at 40 MeV given by the CH89 
OP (dashed line) and mCH89 OP (solid line). The OM prediction for the elastic p+^'^^Pb scattering 
at 54.2 MeV is given by the CH89 OP, with the Coulomb correction to the proton incident energy 
AEc ~ 14.2 MeV. 



with a rather strong AlVc* peaked at the nuclear surface (see Fig. 



IV. DISCUSSION 

We find that p+'^^^Ph data at proton incident energies of 45 and 54.2 MeV have closely the 
same diffraction pattern as elastic neutron scattering at 30.4 and 40 MeV, corresponding to 
energy shifts AEq of 14.6 and 14.2, respectively, owing to Coulomb repulsion effects. Using 
Ep = En + AEc (with the Coulomb correction taken into account by the CH89 systematics) 
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FIG. 9. (Color online) The same as Fig. [7| but for the imaginary part of the mCH89 optical 
potential of 54.2 MeV proton (see Table |II|) . 



we found that the best-fit proton OP still contains a non-zero remnant of the Coulomb 
correction, which represents the isospin impurity of the nucleon OP. While the correction 
to the real part of the proton OP is only a few percent, as might be expected from the 
mean-field nature of the real nucleon OP, the correction to imaginary part at the nuclear 
surface is about 20%. 

Such a significant isospin impurity of the imaginary OP found in our analysis also confirms 
the trend found recently in a global dispersive optical model (DOM) analysis of elastic proton 
and neutron scattering 32]. Namely, from a comparison of the elastic proton and neutron 
data on ^°^Pb, the surface component of the imaginary OP obtained in the DOM has shown 
quite different (A^ — Z) asymmetry dependences between protons and neutrons. Such a 
difference results directly on a deviation of the isospin dependence of the nucleon OP from 
the Lane form (l2])-([3]) and gives rise, therefore, to a larger Coulomb correction to imaginary 
OP at the nuclear surface. 
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These differences indicate that there are different mechanisms for proton and neutron 
absorption that are hnked to different non-elastic reaction channels induced by proton and 
neutron on ^°^Pb target. The most obvious of these are Coulomb excitation and a difference 
in the Q values of the (n, p) and (p, n) reactions. Because the physics origin of the imaginary 
OP is multifaceted and contains dynamic higher-order (beyond the mean-field) contributions, 
the significant isospin impurity found for the imaginary part of the nucleon OP was not 
unexpected. One needs an explicit microscopic calculation of the imaginary OP based, e.g., 
on the Feshbach formalism |31|, to establish full Lane consistency of the mean- field part of 
the nucleon OP. The present results are also in a qualitative agreement with the Brueckner- 
Hartree-Fock calculation of the nucleon OP by Jeukenne et al. jsl, which showed that for 
heavy nuclei such as ^osp]^ ^Yiq imaginary Coulomb correction can be quite significant. 

Our results show that some treatments of the Coulomb correction adopted in the literature 



are probably inadequate. For example, the CH89 systematics [13| uses the same energy shift 
AEc for the real and imaginary OP and it underestimates, therefore, the Coulomb correction 
to the imaginary OP as shown above. The recent global OP by Koning and Delaroche 12 1 
even assumes APFc ~ 0. 

We note that another common usage for Coulomb correction is the entire difference 
of the OP for neutrons and protons at the same energy. To show the strength of such 
a total Coulomb correction, we generated the 45 MeV proton OP using the same CH89 
formulas but setting AE'c = 0. This procedure yields a significantly stronger A?7c, with 
AK/K ~ 5.6%, ^WJW^ ^ 56.9%, and AWs/W^ ^ 12.9%. However, one should be careful 
in discussing such a total Coulomb correction because the CH89 formulas were determined 
[13I using a constant energy shift AEc = 19 MeV for p-t-^^^Pb OP, and it is questionable to 
use the CH89 potential obtained with AEc = in the present discussion. 



V. SUMMARY 

We have measured cross sections for elastic scattering of neutrons from ^osp]^ gg^^ g^j^^ 
40.0 MeV and deduced the realistic OP parameters using the Woods-Saxon geometry given 
by the CH89 systematics. 

The elastic p-|-^°^Pb scattering data at 45 MeV and elastic n-|-^°^Pb scattering data at 
30.4 MeV have about the same diffractive structure at the forward angles, indicating that the 
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energy shift used for 45 MeV protons in the CH89 parametrization should be ^Ec ~ 14.6 
MeV. The isovector part Ui of the nucleon OP for a 208p|^ target has been used in a CC 
analysis of the charge exchange (p, n) reaction at 45 MeV, exciting the IAS in ^°^Bi, and 
a very satisfactory description of the (p, n) data has been obtained. This result allowed us 
to use the complex isovector part t/i of the CH89 OP to check the Lane consistency of the 
nucleon OP. 

The detailed OM analysis of the elastic neutron and proton scattering has shown that the 
realistic proton OP at energy Ep = En + AEc contains a non-zero Coulomb correction to its 
complex strength. Such a non-zero Coulomb correction represents the isospin impurity of 
the CH89 nucleon OP, which is only a few percent for the real part of the OP but is around 
20% for the imaginary part at the nuclear surface. 

We reiterate that these comments are in reference to the CH89 global OP, which already 
contains a significant correction for Coulomb effects. These results show that CH89 system- 
atics provides an accurate estimate of Coulomb effects for the real part of the OP, provided 
that the energy shifts that enter the model are found by matching diffraction structures for 
neutron and proton scattering. The 20% correction for the imaginary part of the OP is not 
large, and the CH89 systematics can still provide a reasonable priori estimate even for the 
imaginary potential. 

These results confirm again the importance of elastic neutron scattering experiments, for 
use in conjunction with existing elastic proton scattering and (p, n) charge exchange data, 
to obtain estimates of Coulomb corrections in heavy nuclei. 



ACKNOWLEDGEMENTS 

We wish to thank Darrell Drake for the loan of the ^osp]^ target. This research was 
supported by the U.S. National Science Foundation grants No. PHY-78-22696, PHY06- 
06007, PHY08-22648(JINA). Two of the authors (D.T.K. and B.M.L.) were supported by 
the National Foundation for Science and Technology Development (NAFOSTED), under 
Project Nr. 103.04-2011.21. 

22 



[1] A.M. Lane, Phys. Rev. Lett. 8, 171 (1962); A.M. Lane, NucL Phys. 12, 676 (1962). 

[2] G.R. Satchler, R.M. Drisko, and R.H. Bassel, Phys. Rev. 136, B637 (1964). 

[3] G.R. Satchler, Direct Nuclear Reactions (Clarendon Press, Oxford, 1983). 

[4] J. Rapaport, V. Kulkarni, and R. W. Finlay, Nucl. Phys. A330,15 (1979); G.R. Satchler, in 

Isospin in Nuclear Physics, ed. D.J. Wilkinson (North-Holland, Amsterdam, 1969), p. 389. 
[5] D.M. Patterson, R.R. Doering, and A. Galonsky, Nucl. Phys. A263, 261 (1976). 
[6] J.P. Jeukenne, A. Lejeune, and C. Mahaux, Phys. Rev. C 15, 10 (1977); 
[7] D.T. Khoa, E. Khan, G. Colo, and N. Van Giai, Nucl. Phys. A706, 61 (2002). 
[8] D.T. Khoa, Phys. Rev. C 68, 011601(R) (2003). 

[9] D.T. Khoa, H.S. Than, and D.C. Cuong, Phys. Rev. C 76, 014603 (2007). 
[10] P.D. Kunz, in The (p,n) Reaction and the Nucleon-Nucleon Force, ed. CD. Goodman, S.M. 

Austin, S.D. Bloom, J. Rapaport, and G.R. Satchler (Plenum Press, New York, 1980). 
[11] F.D. Becheetti and G.W. Greenlees, Phys. Rev. 182, 1190 (1969). 
[12] A.J. Koning and J.P. Delaroche, Nucl. Phys. A713, 231 (2003). 

[13] R.L. Varner, W.J. Thompson, T.L. McAbee, E.J. Ludwig, and T.B. Clegg, Phys. Rep. 201, 
57 (1991). 

[14] E.L. Hjort, F.P. Brady, J.L. Romero, J.R. Drummond, D.S. Sorenson, J.H. Osborne, B. 

McEachern, and L.F. Hansen, Phys. Rev. C 50, 275 (1994). 
[15] J.H. Osborne, F.P. Brady, J.L. Romero, J.L. Ullmann, D.S. Sorenson, A. Ling, N.S.P. King, 

R.C. Haight, J. Rapaport, R.W. Finlay, E. Bauge, J.P. Delaroche, and A. J. Koning, Phys. 

Rev. C 70, 054613 (2004). 
[16] J. Klug et ah, Phys. Rev. C 67, 031601(R) (2003); A. Ohrn et al, Phys. Rev. C 77, 024605 

(2008). 

[17] R.K. Bhowmik, R.R. Doering, L.E. Young, S.M. Austin, A. Galonsky, and S.D. Schery, Nucl. 

lustrum. Meth. 143, 63 (1977). 
[18] R.R DeVito, Ph.D. Thesis, Michigan State University, 1979 (unpubhshed) . 
[19] R.P. DeVito, S.M. Austin, U.E.P. Berg, W. Sterrenburg, and L.E. Young, Nucl. lustrum. 

Meth. 215, 423 (1983). 
[20] Pb target was obtained on loan from D. Drake, Los Alamos. 

23 



[21] W.E. Kinney, Nucl. Instrum. Meth. 83, 15 (1970). 

[22] GIBELUMP was written by F.G. Perey and modified by R.M. Haybron. See also Ref. Q. 
[23] J. Rapaport, T.S. Cheema, D.E. Bainum, R.W. Finlay, and J.D. Carlson, Nucl. Phys. A296, 
95 (1978). 

[24] W.T.H. van Oers, Huang Haw, N.E. Davison, A. Ingemarsson, B. Fagerstrom, and G. Tibell, 

Phys. Rev. C 10, 307 (1974). 
[25] G.S. Mani, D.T. Jones, and D.Jacques, Nucl. Phys. A165, 384 (1971). 

[26] J. Raynal, Computing as a Language of Physics (IAEA, Vienna, 1972) p. 75; J. Raynal, 

coupled-channel code ECIS97 (unpublished). 
[27] R.P. DeVito, S.M. Austin, W. Sterrenburg, and U.E.P. Berg, Phys. Rev. Lett. 47, 628 (1981) 
[28] D.T. Khoa, G.R. Satchler, and W. von Oertzen, Phys. Rev. C 56, 954 (1997). 
[29] R.R. Doering, D.M. Patterson, and A. Galonsky, Phys. Rev. C 12, 378 (1975). 
[30] R.F. Carlson, At. Data and Nucl. Data Tables 63, 93 (1996). 
[31] H. Feshbach, Theoretical Nuclear Physics, Vol. II (Wiley, NY, 1992). 

[32] J.M. Mueller, R.J. Charity, R. Shane, L.G. Sobotka, S.J. Waldecker, W.H. Dickhoff, A.S. 
Crowell, J.H. Esterline, B. Falhn, C.R. Howell, C. Westerfeldt, M. Youngs, B.J. Crowe, and 
R.S. Pedroni, Phys. Rev. C 83, 064605 (2011). 



24 



